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we also show how the shape invariance property of the potential is important in order to 
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some examples of planar vector fields coming from supersymmetric quantum mechanics are 
studied. 
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1 Introduction 

We deal with some generalization of the Darboux theory of integrability for planar vector fields 
and the Darboux transformation of the associated equation. 

In 1882, Darboux in his paper [16] presented as a proposition the notable theorem today 
known as Darboux transformation. This proposition also can be found in his book [17, p. 210]. 
Curiously Darboux's proposition was forgotten for a long time. In 1926 Ince included it in his 
book as an exercise (see Exercises 5, 6 and 7 in [22, p. 132]). Ince follows closely Darboux's 
formulation given in [16, 17]. P. Dirac, in 1930, published The principles of quantum mechanics, 
where he gave a mathematically rigorous formulation of quantum mechanics. In 1938, J. Delsarte 
introduced the notion of transformation (transmutation) operator, today known as intertwining 
operator and is closely related with Darboux transformation and ladder operators. Later on, in 
1941, E. Schrodinger factorized in several ways the hypergeometric equation. It was a byproduct 
of his factorization method originating an approach that can be traced back to Dirac's raising 
and lowering operators for the harmonic oscillator. Ten years later, in 1951, another factoriza- 
tion method was presented by L. Infold and T.E. Hull where they gave the classification of their 
factorizations of linear second order differential equations for eigenvalue problems of wave me- 
chanics. In 1955, M.M. Crum inspired by Liouville's work about Sturm-Liouville systems and 
developed one kind of iterative generalization of Darboux transformation. Crum surprisingly 
did not mention Darboux. In 1971, G.A. Natanzon studied a general form of the transforma- 
tion that converts the hypergeometric equation to the Schrodinger equation writing down the 



*This paper is a contribution to the Special Issue "Geometrical Methods in Mathematical Physics". The full 
collection is available at http://www.emis.de/journals/SIGMA/GMMP2012.html 



2 



P.B. Acosta-Humanez and Ch. Pantazi 



most general solvable potential, potential for which the Schrodinger equation can be reduced to 
hypergeometric or confluent hypergeometric form, concept introduced by himself. 

Almost one hundred years after Darboux's proposition, in 1981, Edward Witten with his 
renowned paper [45] gave birth to the supersymmetric quantum mechanics and he discussed 
general conditions for dynamical supersymmetry breaking. Since Witten's work we can found 
in the literature a big amount of papers related to supersymmetric quantum mechanics. Maybe 
the most relevant of these papers was written in 1983 by L.E. Gendenshtein, where he presented 
the shape invariance condition, i.e. the preservation of the shape under Darboux transforma- 
tion. Gendenshtein used this property to find the complete spectra for a broad class of prob- 
lems including all known exactly solvable problems of quantum mechanics (bound state and 
reflectionless potentials). Today this kind of exactly solvable potentials satisfying the shape in- 
variance condition are called shape invariant potentials, see [20]. In 2009, in [1, 3] were presented 
a Galoisian point of view of the supersymmetric quantum mechanics and in particular of the 
Darboux transformations and shape invariance condition. We point out that the first analysis of 
Darboux transformations from the differential Galois point of view was done by V.P. Spiridonov 
in [40] . The present work follows these approaches with the same point of view. 

From the other hand, Darboux in 1878 presented a simple way to construct first integrals and 
integrating factors for planar polynomial vector fields, see [15]. The key point of his method 
are the invariant algebraic curves of such vector fields. His approach has been related with 
problems concering limit cycles, centers and bifurcation problems, see for instance [21, 31, 38]. 
Moreover, the geometric scenario of the algebraic curves determines the structure of the vector 
fields, see [10, 11, 35]. Nowdays Darboux's method has been improved for polynomial vector 
fields basically taking into account the multiplicity of the invariant algebraic curves see for 
instance [9, 11, 13]. Moreover, the existence of a rational first integral and Darboux's method 
are related by Jouanolou's results, see [23, 32]. Prelle and Singer [36, 39] gave the relation 
between elementary/Liouvillian first integrals and integrating factors that are constructed by 
Darboux's method. Additionally, Darboux's ideas have been extended to a particular class of 
non-autonomous vector fields, see [5, 30]. 

Our main aim in this paper is to relate Darboux's theory of integrability for planar vector 
fields of Riccati type with Schrodinger equation and Darboux transformation of the associated 
equation. So we consider polynomial vector fields of the form 

v = dv/dt = Sq{x) + Si(x)v + S2(x)v 2 , x = dx/dt = N(x), 

with So,Si,S2 E C[x]. In particular for non-relativistic quantum mechanics we have S2 (x) = 
-N(x),Si(x) = and So (a;) = N(x){V(x) - A) with V(x) = T{x)/N{x), N,T G C[x] and A 
a constant. 

Hence, we deal with systems of the form 

v = dv/dt = N(x)(V{x)-X-v 2 ), x = dx/dt = N(x), (1) 
or equivalently we can consider the polynomial vector field 

v = dv/dt = T{x)-N{x)\-N{x)v 2 , x = dx/dt = N(x), (2) 
of degree m = maxjdeg T(x), degN(x) + 2}. Note that the associated foliation of system (1) is 

v' = ^ = V(x)-X-v 2 , (3) 

with V 6 C(x), see also [2]. 

The structure of the paper is the following: In Section 2 we present the basic concepts of 
differential Galois theory, Schrodinger equation, Darboux transformation and Darboux's theory 
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of integrability of planar polynomial vector fields. In Section 3 we present our main results. 
More concrete, in Proposition 3 we characterize the differential Galois group of the Schrodinger 
equation with the class of the first integral of the corresponding vector field. Additionally, in 
Proposition 4 we present a condition for the non-existence of a rational first integral. In Propo- 
sition 5 we show how we can construct the generalized Darboux first integrals and integrating 
factors of the transformed vector field using a solution of the initial Schrodinger equation. In 
Theorem 9 we show that the strong isogaloisian property of the Darboux transformation and 
the shape invariance property of the potential are necessary in order to preserve the rational 
structure of the elements of the transformed vector field. At the end, in Section 4 we give several 
examples as applications of our results. 

As far as we know this is the first time that is presented in the literature the link between Dar- 
boux theory of integrability of planar vector fields and Darboux transformation of the associated 
Schrodinger equation. 



2 Theoretical background 

In this section we present the theoretical background that we use in this work. Most of the 
results are naturally extended in higher dimension. 

2.1 Differential Galois theory 

We start regarding an algebraic model for functions and the corresponding Galois theory known 
as differential Galois theory or also Picard-Vessiot theory, see [26, 27, 37, 43, 44] for all detail. 
The following preliminaries correspond to a quick overview of this theory and can be found also 
in [1, 3]. 



Definition 1 (differential fields). Let F be a commutative field of characteristic zero. A deriva- 

d_ 

dx 



tion of F is a map 4- : F — )• F satisfying 



d , , . da db d , , . da , db 

a + &) = — + —, —a- 6= — -b + a-—, 
dx dx dx dx dx dx 

for all a,b G F. We say that (F, ^) (or just F when there is no ambiguity) is a differential field 
with the derivation 4^. 

We assume that F contains an element x such that 4:{%) = 1. Let C be the field of constants 
of F: C = {c G F | ^ = 0}. C is of characteristic zero and will be assumed to be algebraically 
closed. 

Throughout this paper, the coefficient field for a differential equation will be defined as the 
smallest differential field containing all the coefficients of the equation. 

In particular we deal with second order linear homogeneous differential equations, i.e., equa- 
tions of the form 

d^ d 

^|+a^+/3y = 0, a,(3eF. (4) 

Definition 2 (Picard-Vessiot extension). Consider the differential equation (4). Let L be 
a differential field containing F (a differential extension of F). We say that L is a Picard- 
Vessiot extension of F for the differential equation (4) if there exist two linearly independent 
solutions of the differential equation (4) namely 2/1,2/2 £ L such that L = .F (2/1,2/2) (i- e - L = 
F(yi, 2/2, dyi/dx, dy2/dx)) and moreover L and F have the same field of constants C. 
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In what follows, we work with Picard-Vessiot extensions and the term "solution of (4)" 
will mean "solution of (4) in L". So any solution of the differential equation (4) is a linear 
combination (over C) of yi and 1/2 ■ 

Definition 3 (differential Galois groups). An F-automorphism a of the Picard-Vessiot exten- 
sion L is called a differential automorphism if 

a ( — 1 = —(a (a)) Ma G L and a(a) = a Ma G F. 
\dx J ax 

The group of all differential automorphisms of L over F is called the differential Galois group 
of L over F and is denoted by DGal(L/F). 

Given a G DG&l(L/F), we see that {cyi, 0^2} are also solutions of the equation (4). Hence 
there exists a matrix A a G GL(2, C), such that 

cr{{yi 2/2)) = (cr(yi) cr(y 2 )) = {yi J/2) Ax- 
As <7 commutes with the derivation, this extends naturally to an action on a fundamental solution 
matrix of the companion first order system associated with the equation (4). We have 

( ( yi 2/2 \\ ( °(yi) ^(2/2) \ / yi y2 \ 

a \\dyi dm )) = \ a ( d y±) a( dV2 \) = \ d ^ ^ A<T ' 
\\dx dx I J V \dx) \dx)/ \dx dx/ 

This defines a faithful representation DGal(L/A') — > GL(2,C) and it is possible to consider 
DGal(L/A) as a subgroup of GL(2,C) and depends on the choice of the fundamental system 
{2/1,2/2} only up to conjugacy. 

Recall that an algebraic group G is an algebraic manifold endowed with a group structure. 
Let GL(2, C) denote, as usual, the set of invertible 2x2 matrices with entries in C (and SL(2, C) 
be the set of matrices with determinant equal to 1). A linear algebraic group will be a subgroup 
of GL(2, C) equipped with a structure of algebraic group. One of the fundamental results of the 
Picard-Vessiot theory is the following theorem (see [26, 27]). 

Theorem 1. The differential Galois group DGal(L/-F) is an algebraic subgroup o/GL(2,C). 

In fact, the differential Galois group measures the algebraic relations between the solutions 
(and their derivatives) of the differential equation (4). It is sometimes viewed as the object 
which should tell "what algebra sees of the dynamics of the solutions" . 

In an algebraic group G, the largest connected algebraic subgroup of G containing the identity, 
noted G°, is a normal subgroup of finite index. It is often called the connected component of 
the identity. If G = G° then G is a connected group. 

When G° satisfies some property, we say that G virtually satisfies this property. For example, 
virtually solvability of G means solvability of G° (see [44]). 

Theorem 2 (Lie-Kolchin). Let G C GL(2,C) be a virtually solvable group. Then G° is trian- 
gularizable, i. e. it is conjugate to a subgroup of upper triangular matrices. 

Throughout this work we will use the following definition. 

Definition 4 (Liouvillian integrability) . We say that the linear differential equation (4) is 
(Liouville) integrable if the Picard-Vessiot extension L D F is obtained as a tower of differential 
fields F = Lq C L\ C • • • C L m = L such that Li = Lj_i(r/) for % = 1, . . . , m, where either 

1) 7] is algebraic over that is r\ satisfies a polynomial equation with coefficients in 
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2) 7] is primitive over that is G 



3) ?7 is exponential over that is 




We remark that in the usual terminology of differential algebra for integrable equations the 
corresponding Picard-Vessiot extensions are called Liouvillian. From now on we say that an 
equation is integrable whether it is integrable in the sense of differential Galois theory according 
to Definition 4. The following theorem is due to Kolchin. 

Theorem 3. The equation (4) is integrable if and only if~DG&l(L/F) is virtually solvable. 

There is an algorithm due to Kovacic [28] that decides about the integrability of the equa- 
tion (4) in the case where F = C(x). In practice, Kovacic's algorithm deals with the reduced 
form of equation (4), namely with the form y" = ry, where r is a rational function. Kovacic 
used the fact that DGal(L/F) C SL(2,C), in order to separate his algorithm in three cases for 
the integrability of the equation (4): 

Case 1. DGal(L/F) is reducible, 
Case 2. DGal(L/F) is irreducible, 
Case 3. DGal(L/F) is finite primitive. 



2.2 The Schrodinger equation 

Here we first introduce the Schrodinger equation and then we present the preliminaries about 
Schrodinger equation from a Galoisian point of view, see [1, 3]. 

In classical mechanics for a particle of mass m moving under the action of a potential U the 
Hamiltonian is given by 

IIpII 2 

H = — VU(x), p = (p 1 ,...,p n ), x = {xi,...,x n ), 

2m 

and corresponds to the energy (kinetic plus potential). From the other hand in quantum me- 
chanics the momentum p is given by p = —ihV, where h is the Planck constant and V is the 
Laplacian operator. In this case the Hamiltonian operator H is the Schrodinger (non-relativistic, 
stationary) operator which is given by 

H = -fv* + V(x), 
2m 

where x is the coordinate and V (x) is the potential or potential energy. The Schrodinger equation 
is given by = A^, where the eigenfunction \t is the wave function and the eigenvalue 
A € Spec(-ff) is the energy level. The solutions $ of the Schrodinger equation are the states of 
the particle and Spec(-ff) denotes the spectrum of the operator H. In [41] it can be found the 
details about the mathematical foundations of quantum mechanics for the Schrodinger equation. 

According to [14, 45], a supersymmetric quantum mechanical system is one in which there 
are operators Qi that commute with the Hamiltonian % and satisfying 



[Qi,H] = QiH - HQi = 0, 
xQij Qjj — QiQj QjQi — ^%0~l 



and 



1, i = J, 
0, i^j. 
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For n = 2, we obtain the simplest example of a supersymmetric quantum mechanical system. 
In this case we have that Thus, the supercharges Qi are defined as 

_ a 1 p±a 2 W(x) 
Q± — 2 ' V+ — Qii Q-—Q2, 

where p = W : R — > R is the superpotential and <7j are the Pauli spin matrices. 

The operator satisfying = HQi and 2Q? = H, is given by 

I 2 p 2 + I 2 W 2 (x) + ha 3 £W(x) (H + \ _ /l 

H = 2 = V0 ffj' l2= V0 1 

The operators ff_ and i7 + are the supersymmetric partner Hamiltonians and are given by 

where V± are the supersymmetric partner potentials. 

Now we follow [1, 3] in order to present the Schrodinger equation in the context of differential 
Galois theory. Thus, the Schrodinger equation (stationary and one dimensional) now is written as 

d 2 

HV = \y, H = + V(x), VeF, (6) 

where F is a differential field (with C = C as field of constants). We will deal with the 
integrability of equation (6) in agreement with our definition of integrability, i.e., in the sense 
of differential Galois theory, see Definition 4. In [1, 3] were introduced the following notations, 
useful for our purposes. 

• A C C denotes the algebraic spectrum of H, i.e., the set of eigenvalues A such that equa- 
tion (6) is integrable (Definition 4). 

• L\ denotes the Picard-Vessiot extension of equation (6). Thus, the differential Galois 
group of (6) is denoted by DGal(LA/-^)- 

Definition 5 (algebraically solvable and quasi-solvable potentials). The potential V(x) G F is: 

• an algebraically solvable potential when A is an infinite set, or 

• an algebraically quasi- solvable potential when A is a non-empty finite set, or 

• an algebraically non-solvable potential when A = 0. 

When Card(A) = 1, we say that V{x) £ F is a trivial algebraically quasi-solvable potential. 

The following theorem shows that if there exist more than one eigenvalue in the algebraic 
spectrum of the Schrodinger operator with F = C(x), then we cannot fall in case 3 of Kovacic's 
algorithm. 

Theorem 4 (see [1, 3]). Consider the Schrodinger equation (6) with F = C(x) and Picard- 
Vessiot extension L\. 7/DGal(Lo/-F) is finite primitive, thenT)G&l(L\/F) is not finite primitive 
for all A G A\{0}. 

From [1, 3] note that the known cases of rational potentials in quantum mechanics leads to 
Schrodinger equations falling in case 1 of Kovacic's algorithm. Additionally, if Card(A) > 1 
then any algebraic solution of the Riccati equation associated to the Schrodinger equation (6) 
is a root of a polynomial of degree at most two. 
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2.3 Darboux transformation 

Darboux gave in [16] a transformation that allow us to transform some type of differential equa- 
tions into other differential equations preserving the type. The following results corresponds to 
the Darboux transformation, denoted as DT, in the Galoisian and quantum mechanic formalism, 
see [1, 3]. 

Theorem 5 (Galoisian version of DT). Assume H± = —-^z + V±(x) and A / 0. Consider 
the Schrodinger equation H-^~^ = X^Sf^ with V—(x) G F. Let DT be the transformation such 
that V- h-> V + , VlK - ) (->■ \I>( + ), Fh-F. Then for the Schrodinger equation H + ^ +s > = A\I/( + ) with 
V-\-{x) G F the following statements holds: 

DT(V_) _ V + _ + A, _ V_ - 2^(ln*H), 

DT(^ ( ~ } ) = ^ +) = — w/iere 

^ is a particular solution of = X\^^\ Ai G A; 

ii) DT(¥<->) = *W = iU<-> _ A(i n *H)*H A ^ Alj w/iere 

\|A is i/ie general solution of H-^^ = A 1 ^"- 1 /or A G A \ {Ai} and 
^ +) is i/ie general solution of H + ^ + ^ = also for A G A \ {Ai}. 

Remark 1. According to [1, 3] we have that a transformation is called isogaloisian whether it 
preserves the differential Galois group: the initial equation and the transformed equation have 
the same differential Galois group. Furthermore, when the differential field and the Picard- 
Vessiot extension are preserved, then the transformation is called strong isogaloisian. 

In agreement with Theorem 5 and Remark 1, we obtain the following results, see [1, 3] for 
complete statements and proofs. 

Proposition 1. In general, DT is isogaloisian and virtually strong isogaloisian. Furthermore, 
if dxO- n ^\-^) S F, then DT is strong isogaloisian. 

Proposition 2. The supersymmetric partner potentials V± are rational functions if and only if 
the superpotential W is a rational function. 

Corollary 1. The superpotential W G C(x) if and only i/DT is strong isogaloisian. 

Remark 2. The examples of this paper are in agreement with the previous results, since F = F 
due to the fact that the superpotential W belongs to F = C(x). The following definition is 
a partial Galoisian adaptation of the original definition given in [20] (F = C(x)). The complete 
Galoisian adaptation is given when F is any differential field, see [1, 3]. 

Definition 6 (rational shape invariant potentials, see [1, 3]). Assume V±(x; fi) G C(x; /u), where 
[i is a family of parameters. The potential V = V- G C(x) is said to be rational shape invariant 
potential with respect to \i and A = A n being n G Z + , if there exists a function / such that 

n+l 

V+(x;a ) = V-(x;ai) + R(a{), a 1 = f(a ), A n = ^i?(a fe ), A = 0. 

fc=2 

Hence the form of the potentials V± are preserved up to parameters. Theorem 5 and Propo- 
sitions 1 and 2 lead us to following result. 
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Theorem 6 (see [1, 3]). Consider Hty( ' = \ n ^( ) with Picard-Vessiot extension L n , where 
n G Z+. IfV = V- £ C(x) is a shape invariant potential with respect to A = X n , then 

DGal(L n+ i/C(x)) = DGal(L n /C(aO), n > 0. 

Hence for rational shape invariant potentials Galois group is preserved due to Picard-Vessiot 
extension and differential field. 



2.4 Darboux's theory of integrability for planar polynomial vector fields 

In this subsection we present the basic ideas of Darboux's method for planar polynomial vector 
fields, see [15]. We don't give en extensive presentation of this theory but we only present the 
basic results that we need in Section 3. 

We consider the polynomial {differential) system in C 2 defined by 

~ = i = P(x,y), ~=y = Q{x,y), (7) 

where P and Q are polynomials in the variables x and y. The independent variable t can be real 
or complex. We associate to the polynomial differential system (7) the polynomial vector field 

x = P(x,y)^ + Q(x,y)^, (8) 

and its associated foliation is given by Pdy — Qdx = 0. 

An algebraic curve f(x,y) = in C 2 with / 6 C[x,y] is an invariant algebraic curve of the 
vector field (8) if 

A(log(/)) = K, (9) 

for some polynomial K £ C[x, y] called the cof actor of the invariant algebraic curve / = 0. Note 
that due to the definition (9) we have that the degree of the cofactor K is less than the degree 
of the polynomial vector field (7). Moreover, the curve / = is formed by trajectories of the 
vector field X. 

For a given system (7) of degree m the computation of all the invariant algebraic curves 
is a very hard problem because in general we don't know about the maximum degree of such 
curves. However, imposing additionally conditions either for the structure of the system or for 
the nature of the curves we can have an evidence of a such a bound [7, 8, 12]. 

Let h,g G C[x, y] be relatively prime in the ring C[x,y]. The function exp(g/h) is called 
an exponential factor of the polynomial system (7) if there is a polynomial L £ C[x,y] (also 
called cofactor) that satisfies the equation X (g/h) = L. It turns out that if h is not a constant 
polynomial, then h = is one of the invariant algebraic curve of (7). 

The following theorem is a short version of Darboux theory of integrability for planar poly- 
nomial differential systems. For more details and also for higher dimension see [29]. For genera- 
lizations see [5]. 

Theorem 7. Suppose that a polynomial system (7) of degree m admits 

• p irreducible invariant algebraic curves fi = with cof actors Ki for i = 1, . . . ,p, 

• q exponential factors Fj with cof actors Lj for j = 1, . . . , q. 



Then the following statements hold: 
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(a) The function 

f^---fp p F^---F^, (10) 

p <? 

is a first integral of (7) if and only if ^ \%Ki + ^ ^jLj = 0. 

1=1 i=l 

(6) T/ie function (10) is a (Darboux) integrating factor of the vector field (8) i/ and on/y z/ 

jr \iKi + s + (p x + Q y ) = o. 

i=l i=l 

Integrating factors for planar systems can be thought as parametrization of the independent 
variable (time) that yields to a divergence free system. 

As we show in Theorem 7 for polynomial vector fields the invariant algebraic curves and 
(because of their multiplicity, see [13]) the exponential factors are the basic elements in order 
to construct first integrals/integrating factors. As we will see the method of Theorem 7 also 
works for more general expressions of vector fields, curves, exponential factors and cofactors, 
see also [5, 18, 19]. 

Definition 7 (generalized exponential factor). Consider S(x) £ C(x). We define a generalized 
exponential factor of a polynomial vector field X any expression of the form F = exp (J S(x)) 
which satisfies X(F) = LF and L is called generalized cof actor. 

Definition 8 (generalized Darboux function). We define generalized Darboux function any 
expression of the form 



(y - St(x))^ • ■ ■ (y - S p (x)yv exp U S(x) 

= (y- S 1 (x)) x ^ ■■■(y- S p (x)) x ? exp (S(x)) JJ(x - (11) 

with Si, S,S,g S C(x), Xi, hi, Aj € C and p G Z+, see also [46]. 

An integrating factor of the form (11) will be called generalized Darboux integrating factor 
and similarly a first integral of the form (11) will be called generalized Darboux first integral. 

Remark 3. We consider the Schrodinger equation (6) with F = C(x) and the associated 
Ricatti equation (3). We also consider v\(x), V2(x), i>3(x) particular solutions of equation (3) 
with V £ C(x) and we write V(x) = T(x)/N(x) with T,N £ C[x]. Note that an associated 
polynomial vector field to equation (3) can be written into the form (2). 

(a) Following [2, 28, 46] we can distinguish the following cases about the type of the first 
integrals of equation (3) or equivalently of the polynomial vector field (2). 

Case 1. (i) If only v\(x) £ C(x) then the vector field (2) has a first integral of Darboux- 
Schwarz-Christoffel type, namely a first integral of the form 

X 

I(v, x) = exp(o(x)) TT(x - xi) a * + / exp{g(u)) TT(u - Xi) ai ~ mi P{u)du, 

y-vi(x) J 11 

with P £ C[x], g £ C(x), Xi,a{ £ C and rrii £ Z + , see also [46]. 

(ii) If both Vi(x), v%(x) £ C(x) then the vector field (2) has a generalized Darboux first 
integral of the form 

I(v,x) = ~ v + v ^\ exp ( [ ( V2 ( x )- Vl ( x ))dx) . (12) 
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Furthermore if the differential Galois group of the Schrodinger equation (6) is a cyclic group of 
order k then the k-ih. power of the first integral (12) is rational, see also [2]. 

Case 2. If v% is a solution of a quadratic polynomial then the vector field (2) has a first 
integral of hyperelliptic type. 

Case 3. If all v\, V2, V3 are algebraic over C(x) then the vector field (2) has a rational first 
integral of the form 

. , (V2 ~ Vl)(vi ~ V) 

I{V,X) = r. 

(V 3 - vi){v 2 ~ V) 

(6) In general, knowing one algebraic solution v\{x) of equation (3) we can obtain a second 
solution of (3), namely 

^ exp (—2 J v\dx) 
f exp (—2 f v\dx) dx 

Then, the vector field (2) has always a first integral of the form (12) and can be rewritten either 
as Darboux-Schwarz-Christoffel type or as a generalized Darboux function or as a first integral 
of hyperelliptic type or as a rational first integral. 

Additionally in [2] appears the following result about the existence of a rational first integral 
of a polynomial vector field whose foliation is of Riccati type. 

Theorem 8. The Galois group of the equation y" = r(x)y with r G C(x) is finite if and only if 
there exist a rational first integral for the associated polynomial vector field of the corresponding 
Riccati equation. 

Note that the Schrodinger equation (6) can be always written into the form y" = r(x)y and 
so we can always apply Theorem 8. 

3 Main results 

From now on we consider the Schrodinger equation (6) with potential V = V-\ 
*(-)" = (V-(x) - X)^ { -\ V - = J[> T,NeC[x], 

A G A, Card(A) > 1, F = C(x). (13) 

After the change of coordinates ) = ), equation (6) can be written as 

((-> = y_ — \ - C (-)a > (14) 

and we associate to equation (14) the polynomial vector field 

xP-iT-W-N^-^ + N^L, (15) 

of degree m = max{degT(x), deg N(x) + 2}. By Propositions 1, 2 and Corollary 1, we have 

that for given V- G F and W = Co } = ^0 ^/^o } e F then we obtain that V+ G F, i.e., the 
Darboux transformation DT is strong isogaloisian. Moreover, the superpotential W is rational. 
The applications considered in this work correspond to this case. 
The following proposition follows directly by Remark 3 and [2] . 
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Proposition 3. The differential Galois group DGal(L^/ F) of the Schrodinger equation (13) is 
virtually solvable if and only if the first integral of the vector field (15) can be written in one of 
the forms appearing in Remark 3(a). 

In the next proposition we present a result about the non-existence of a rational first integral. 

Proposition 4. Consider the Schrodinger equation (13). IfT)G&l(L\/F) is not cyclic then the 
associated foliation of the Schrodinger equation (13) has not rational first integrals. 

Proof. It follows directly from Theorems 4 and 8. ■ 

Lemma 1. Consider \x) a solution of the Schrodinger equation (13) and we denote by 




(16) 



Then for all X £ A the vector field (15) admits the following: 



• Invariant curve 



/ (-)( C H,x)=-cH+cr ) (x), 



with generalized cofactor 



• Generalized exponential factor 




with generalized cofactor 



4 ) {((-\x)=N'(x)/2 + N(x)(l \x). 



• An integrating factor of the form 





• A first integral of the form 




with C(A,i) as i n (16) and 



■(-) 

(A,2) 



■(-) 
(A,l) 



+ 
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Proof. Note that ({ \x) is a solution of the associated equation ) = V-(x) — A — ^ for 
all A € A. Using the expression of the vector field (15) we have 

xV{fx\^\ X )) = (T-\N-N(^)(-1)+N(^' 

Note that the curve f[ is polynomial in the variable C^ - -* but it could be not polynomial in 
the variable x. Direct computations shows that system (14) admits the generalized exponential 
factor 

FM(cH, I )=e X p(/(i^M + <W(x))<ix), 

with generalized cofactor \^~\x) = N'(x)/2 + G \x). Note that vector field (15) has 
divergence div(X A ) = -2iVC (-) + N'(x). Hence, we have that -21^ ) - 2L^ } + div(X^ } ) = 
and similarly to Theorem 7(6) we have that i? A ^ = (f^ ^ f[ h 2 is an integrating factor for 
the vector field (15). 

If ty^^Jx) a solution of the Schrodinger equation (13) then K is a solution of the Riccati 

equation (14). Then, according to Remark (3)(6) we have that L is another solution of the 
Riccati equation (14). In this case the vector field (15) admits the two invariant curves 

&1) (c H ,*) = + 41)^' = - c(_) + 

with generalized cofactors 

*&)(C H ,*) = -iV(x)(cH + C f- 1) (x)), ^(C^.x) = -JV(x)(cH + C ( ( A 1)(x)). 

The two generalized exponential factors 

FH, ((<->,,) =exp(/ Q^M tjH,,))^ 

with generalized cofactors 

Note that it holds — Kk L + lO^ L ~~ -41) + ^(A2) = ^ anc ^ s i mnar to Theorem 7(a) we have 
that the vector field x{ admits the first integral 

A-) p(-) 

,(-) _ J(\,2) • r (A,2) 

W f (-) p(-) • 

•W) ' *(A,1) 

Hence, the lemma is proved. ■ 

Proposition 5. Consider \x) a particular solution of the Schrodinger equation H- 1 ^!^ = 
\yfr(~) for A = and let (q \x) = (in (\I/q \xj)Y . Then after the Darboux transformation DT 
for all A 7^ equation (14) becomes 

((+)' = V+ - A - C {+)2 , (17) 
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with V+ = — Co^ ^' + Co*" ^ and its associated vector field is 

4 +) = (-d- ) '+Ci- )2 -A-CW 2 )^ + | ; . (18) 

The vector field (18) /or all X ^ admits the following: 

• Invariant curve 

A +) (c (+ \x) = - c « + ci- } + (in(c^ -cS _) ))', 

■urei/i generalized cofactor 

4 +) (C (+) ^) = -C (+) -cl- ) -(ln(cl- ) -C^ ) )) / . 

• Generalized exponential factor 

^ +) (c (+ ^) = (cr ) -ci- ) )exp(/ C |- ) ), 

■urei/i generalized cofactor 

4 +) = <i- ) + (i-(d- ) -rf- ) )) / - 

• ^4n integrating factor (for A 7^ 0) 



i?i +) (C«,x) 



exp ( -2 / Ci _) 



^4 /irsi integral of the form 



i^)-(/(ca-ca)^ 



Proof. For A = we consider \&q ^ a particular solution of the Schrodinger equation H-ty( ) 
A\I/(~). We denote by Co = (l n (^0 ))'■ Then from Theorem 5 for A / we have 



*«=*<->'- <<->*<->= exp (/ ( r>) (cr»-d-'), 



and we have use that (J; ) = *i } - Note that 

is a solution of the Schrodinger equation 

= A^(+). Let 



; In (in (exp (/C^) (C^ - C^))' • 



According to Theorem 5 we have 

^(-^)W>(-§)W> + <r, 
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and we additionally consider the change of variables = (ln^W))'. Then the Schrodinger 
equation H+^S?^ = A\lK + ' falls in equation (17). 

Therefore, similar to Lemma 1, the vector field (18) for all A ^ admits the invariant 
curve f^ with generalized cofactor If j: given by 

/l +) (c (+ U) = -c (+) + = -c {+) +c A - } + (in(cl- } - c^))', 
4 +) (c (+) ,x) = -c« - ci +) = -c (+) - ci- } - (in (d" } - 

Additionally, for all A 7^ admits the generalized exponential factor 



W(C«,x) = exp (/ ci +) ) = exp (/ (C« + (In (C^ - C ^))'] 

with generalized cofactor = = C A ~ } + ( In (Cjf 5 - Co _) ))'- Note that the vector field (18) 
has divergence divA = -2C (+) and for all A ^ it holds -2A^ +) - 2L^ +) + divA A (+) = 0, and 
similarly to Theorem 7(6) we have that = f\^F\( + ^ 

Note that if f is a solution of the Schrodinger equation H+ty^ = A^^ then Cn["i) is 

a solution of the Riccati equation C\ + ^ = V+ — A — £( + ) 2 . Following the same arguments as 
in the proof of Lemma 1 we obtain the expression of the first integral. Hence, the proof is 
completed. ■ 

Theorem 9. If Darboux transformation is strong isogaloisian for a potential V- £ C(x), then 
each pair of functions f^iC^^x), A^^W, 2;), L^\^ ±s, ,x) and I^\c^\x) belong to the 
differential extension in the variables and x. Furthermore, ifV- is a rational shape invariant 
potential, then /j^ (C^> x) , K^(J^\ x) and L^\^^\ x) are rational functions in andx. 

Proof. The proof follows directly by Theorem 6 and Proposition 5. ■ 



4 Applications in quantum mechanics 

In this section we present some examples that relate polynomial vector fields, according to 
Lemma 1 and Proposition 5, with supersymmetric quantum mechanics. Schrodinger equations 
with potentials such as free particle (explicitly discussed in [40] with a similar approach), three 
dimensional harmonic oscillator and Coulomb potential are analyzed as Schrodinger vector fields. 

Example 1 (free particle). Consider the Schrodinger equation H_\l/(~) = A 1 ^ - - 1 with potential 
V- = and differential field 1 F = C(x). Thus, we have that A = C. 

If we choose Ao = and as particular solution of the Schrodinger equation the solution 
= X we have that ci _) = *f," ) 7*o" ) = V^- Since Cq £ C(x) from Proposition 1 we have 
that DT is strong isogaloisian. Additionally, 

and so V- is not shape invariant. Then for A / the general solution of H—^~^ = A\l/( - ) is 
given by 

^ ) =C^(A, ) l)+C2^, ) 2 ), 

x Although the smallest differential field containing the coefficients of the Schrodinger equation is C, to avoid 
triviality, we choose C(x) as the suitable differential field for such equation. 
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>(-) 



where we have consider ^ 



(-) 



,^ — exp [V— Ax) an d ^(A 2) = ex P ( — ^ ~ ' Hence we have that 
ci exp [y— Ax) — C2 exp (— \J— Ax) 



Additionally, we have that 



ci exp (V— Ax) + C2 exp (— \J — Ax) ' 



-A 



and 



(-) 

(A,2) 



-A. 



Here in the expression of the vector field (15) we have that T(x) = and iV(x) = 1. Hence 
the polynomial vector field (15) associated to the Schrodinger equation H—ty^' = X^' with 
potential V- = 0, is 

and is quadratic. Then according to Lemma 1 for all A G A \ {0}, the vector field x| ^ admits 
the following: 
Invariant curve 



ft ) (c i ' ) ^) = -c { - ) + 



(-) 



-ci exp (v 7 — Ax) — C2 exp (— y 7 — Ax) 
ci exp (V— Ax) + C2 exp (— v 7 — Ax) ' 



and f[ ((( \x)gC(C ( - ; , x, exp (v 7 — Ax)) with generalized cofactor 

ci exp (V— Ax) + C2 exp (— \J — Ax) 

and x) G C(C ( ~ } , x,exp (x/^Ax)). 

Generalized exponential factor 



\ x) = c\ exp (v 7 — Ax) + C2 exp (— \J — Ax) , 
with generalized cofactor 



c\ exp (V— Ax) + C2 exp (— \J — Ax) 

Generalized Darboux integrating factor it!^ '(£(~),x) given by 

1 



G C(C (_) ,x,exp (V^Ax)). 



(ci exp (a/^Ax) + c 2 exp (--y/^Ax)) |— ) + \ 

and note that R^\c ( ~\x) G C(C (_) ,x,exp (V 11 ^)). 
Generalized Darboux first integral 



-Ci exp (V— Ax) — C2 exp (— v 7 — Ax) 
ci exp (v 7 — Ax) + C2 exp (— v 7 — Ax) 



2 > 



: exp 



-CH+V-A 
and 4 _) (C (_) ,x) G C(C ( " ) ,x,exp( v /z Ax)). 



-Ax). 



A / 
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Applying the Darboux transformation DT for A ^ (and according to Proposition 5) we 
have 



DT(\^ -) ) = * 



(+) _ ^(~)' _ A-)m(-) 



ci(y/—Xx — 1) exp (y/— Ax) ci{\[— Ax + 1) exp 



-Ax) 



,(+) 

'A 



with 



A = c\ exp (v 7 — Ax) + C2 exp (— v 7 — Ax) 
Additionally we have 



XAx 2 + J^XBx - A 
x{V^XBx - A) ' 

B = c\ exp (V— Ax) — C2 exp (— v— Ax) . 



A+) 



(+) 

(A,2) 



-A + I In 
/ =A + | In 



-A + 



x 



-x- 1 - 

X 



1 

-Ax 2 — x ' 
1 



-Ax 2 + 



We can see that for all A G A the Picard-Vessiot extensions are given by Lq = Lq = C(x), L\ = 
L\ = C (x, exp(\/— Ax)) for A G C*. In this way, we have that DGal(Lo/^) = DGal(Lo/-F) = e; 
for A / 0, we have DGal(L A /F) = DGal(L A /F) = G m , see [1, 3]. 

Now according to Proposition 5 we can compute DT(X A ~) = and we obtain the rational 
vector field 



A 



d d 
+ 



d£(+) dx 



VA ^ 0, 



and equivalently we can consider the polynomial vector field of degree four 



AN 



2 - Ax^ - x z C 



d 



+ x z 



d_ 

dx ' 



VA / 0, 



and admits the following: 
Invariant curve 



7i +) (C (+) ^) = /| +) (C (+) ,x) = - C w + cl- ) 



+ 



A-)', i 
_ I 

'A a; 



-Ax 



cit4 2 exp (a/^Ax) + c 2 i?2 exp (- 
(c\Ai exp (v 7- Ax) — C2-B1 exp (— V— Ax)) ' 



with 



A 1 = V^Xx-l, B 1 = V^Xx + l, A 2 = Ax 2 + Ai, B 2 = Ax 2 - B 1 , 
and note that f^\(,( + \x) € C((( + \ x, exp(\/— Ax)) with generalized cofactor 

Ai +) (C (+) ,x) =x 2 a1 + >(c« x) =x 2 f-C (+) -ci 



<-(-) ^A + ^ 



''A x 



C (+) + 



ciA 2 exp (V— Ax) + C2-B2 exp (— \J — Ax) 
(ciAi exp (v 7 — Ax) — C2-B1 exp (— V— Ax)) 



and note that K { x +) £ C(C (+) , x, exp( v / -Ax)) . 
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Generalized exponential factor 

F[ + \C^,x) = Fi + \c i+ \x) = (C« - exp (/ C. 

c\A\ exp [y— Ax) + C2B2 exp (— \J — Ax) 



(-) 

A 



with generalized cofactor 



(+) 



2r(+) _ „2^(+) 



x z L 



(-) , + g 



X-) _ I 

x 



3 ci^4 2 exp (v^Ax) + c 2 B 2 exp (— V— A) 
ci^4i exp (V— Ax) + C2-B2 exp (— \J — Ax) ' 

and note that € C(C*- + - 1 , x, exp(\/— Ax)) . 
Generalized Darboux integrating factor 



4 +) (c (+) , 



i?i +) (C«,x) 



(ci exp (V— Ax) + C2 exp 



-Ax))" 



- 2 (c^ - j 



-c (+) + c[- ] + 

x 



and E^ +) (C (+) ,x) G C(C (+) ,x,exp(V^Ax)). 
Generalized Darboux first integral 



-c (+) 



1 

1B1 



A /Bi 



exp (— 2V— Ax) , 



and 7j +) (C (+) ,x) G C(C (+) , x,exp(y /= Ax)) 

»(-) 



We could consider £^ L 



obtain the Schrodinger polynomial vector field 
x { ~ ] - (-X - C^ 2 )- 

^(A,l) - { A ^ ) 



AtanVAx, being ^/^ i) = cosVXx and ^ = x. Thus we 



d d_ 

'dC^y + ax' 



of degree two. According to Lemma 1 for all A G A \ {0}, the vector field ^ admits the 
following: 

Invariant curve 

/fci)(C H '^) = -C (_) - v / Atan(v / Ax) G C(C {_) ,x,tan (VAx)), 
with generalized cofactor 

K[~\j((-\x) = -C (_) + v / Atan(v / Ax) G C(C ( ~\ x, tan (VAx)) . 
Generalized exponential factor 



F { {x ^(C ( ~\x) =exp Jtan(y/X. 



Xx)dx 



1 + tan 2 (\/Ax) 



with generalized cofactor ^(C^ \ x ) = — \/Atan(\/Ax) G C(£( ), x, tan(VAx)) . 
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Generalized Darboux integrating factor (for A ^ 0) 

tfcU<*-\x) = 1 + tan2 (^) e C^Utan^)). 

' C (_) - \/Atan(\/Ax) 



Generalized Darboux first integral 

41)(C ( -U) = G C( C « Man^x)). 



(A ' 1)V ' ; (CH + VXtan(v^A)x))VA 
By Proposition 5, Darboux transformation of the vector field XZ ^ is given by 

and we can work with the polynomial vector field of degree four 
X + - (2- \x 2 - x 2 C^ 2 ) 9 +x 2 — 

A (A,1) - [ Z AX X <> )q^+) + X Q x i 

that admits: 

The invariant curve 

f+ ^ = -ti+) 4. -^ x ^(^f\x)+\x 2 -I 

J ^ ,} ~ ^ x(VAa:tau(VAa:)+l) ' 

with generalized cofactor 

-\/Ax tan (\/Ax) + Ax 2 — 1 \ 



^(CW *)=* 2 -c (+) 



.T 



The generalized exponential factor 



(\/Axtan (VAx) + l) 



- . / 1 \ — vAtan (vAx) — - 

^a,i)^ ' x >- r V77T \ ' 

J l + tan 2 (vAx) 

with generalized cofactor 

_ + . . — \/Ax tan (\/Ax) + Ax 2 — 1 
L ^ {C ,X)=X v / Axtan(VAx) + l " 

The integrating factor 

= 1 + tan ' (v/Il 2 ) 

^(-c« + ^) 2 (-f) 2 

with ^3 = \/Ax tan(\/Ax) + 1. 

Also we could consider the partial solution ^r\2) = srn (V^\) and we could work in a similarly 
way as before. 
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Remark 4. As in the same philosophy of the original Darboux transformation, we can iterate 
it to obtain families of new potentials. Starting with V = 0, the following potentials can be 
obtained using Darboux iteration DT n (see [4, 6, 22, 24, 25, 33, 42]): 

_n(n-l)b 2 . m 2 n(n - l)(b 2 - a 2 ) 

1 ) *n — ~77 , ) J *n 



(bx + c) 2 ' n (acosh(ma;) + frsinh(mx)) 2 

— Aabm 2 n(n — 1) m 2 n(n — l)(b 2 + a 2 ) 



(aexp(mx) + 6exp(— mx)) 2 ' ' n (acos(mx) + 6sin(mx)) 2 

In particular for the rational potential given in I), see [1, 3], we have F = F n = C(x) and for 
A„, = A = 0, we have 

= jj-^- + c 2 (bx + cT+\ 
(bx + c) n 

so that 

DGal(L /F) = DGal(4 n) /F) = e, 
whilst for A ^ and A n = 0, the general solution vl/^ is given by 

$W ( x ) = ci A n (x, A)C n ( sin (VAx) + c 2J B n (a;, A)L>„ ( cos (VXr) ) ) , 

where A n , B n , C n , D n £ C(x), so that 

DGal(L A /F) = DGal(4 n) /^) = G m . 

Finally, considering 6 = 1, n = £ + 1, c = and i = rwe arrive to the square well potential, 
which can be considered in a similar way. 

Example 2 (3D Harmonic oscillator). We take V- = r 2 + — {2£ + 3), following [1, 3] we 

see that A n Spec(-ff) = 4N. For instance we consider A = An, where n £ N. Additionally, we 
have 

= r^ 1 exp (-£) , rf"> = ^ - r G C(r), = r^P n exp 

,(-) = t±A + P k_ r A-) _ exp(r 2 ) / exp(r 2 ) j* £ + 1 

\n,l) r P n ' S («. 2 ) r 2£+2p2 y r 2M-2p2 " x " r ^ r ' 

where P n are polynomials of degree n which are related to the generalized Laguerre polynomials, 
see [34]. According to Proposition 1 we see that DT is strong isogaloisian. Since, 

DT(y_) = V + = r 2 + {( + *W + 2) - (2£ + 1), 

we have that the potential V- is also shape invariant. 

Note that T(r) = r A + £(£ + 1) — (2£ + 3)r 2 and N(r) = r 2 . Hence, system (15) can be written 

as 

X- = (r 4 + £{£ + 1) - (2£ + 3)r 2 - 4nr 2 - r 2 Q^ 2 ) -JL- )+r 2 ^-. 
According to Lemma 1 admits the invariant curve 

ftK&Kr) = -C^ + — + #"-£ C(C ( -\r), 
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with generalized cofactor 

K { n ] {C { -\r) = -r 2 C H -{£+ l)r - §r + r 2 G C(C^\r) . 

Ml 



The generalized exponential factor 



with generalized cofactor 

pi 

L { n>{( { ~\r) = r+(£ + l)r + -^r 2 -r 3 G <C(r). 

Ml 

Moreover, admits the generalized Darboux integrating factor 

exp (r 2 ) 



4-)(C ( - } ,r) 



p2 r 2«+2) _ C (-) + |±1 + 



1 , ^ 



r 



where it^ G C(r, exp (r 2 /2)) and the first integral 

j(-)( C (-) r ) = exp ( r2 ) + / exp ( r2 ) dr e j? 

with A\ = r 2 + £( _ V — £ — 1. According to Proposition 5 using the Darboux transformation we 
obtain the rational vector field 

x ^ (r2+ ™ 2 )_ (2f+1) _ 4 „_ <(+)2) ^ + , i 

or equivalently we can consider the polynomial vector field 

X+ = (r A + (£+l)(£ + 2)- (21 + l)r 2 - 4nr 2 - r 2 C (+)2 ) ^ + r 2 |-. 
Moreover, we have 

A+) - P n - | 4n + 1 
\n,l) - pi r + r ' 

-•■1-1 ' 



C = r^ a P„« - (r 2 - I - l)P' n ) J + exp (r 2 ) (r 2 -l-l). 



The vector field admits the invariant curve 

fi +) (C (+ \r) = -C W + ^ + ^ - r G C (C« , r) , 



with generalized cofactor 

(C (+) , r) = -r 2 C (+) - r(* + 1) - ^r 2 + r 3 G C(C (+) , r) . 



The generalized exponential factor 



P' 
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with generalized cofactor 

4 +) {( {+ \r) = r + (£+l)r- ^r 2 + r 3 G C(r). 

n 

Additionally, the vector field admits the generalized Darboux integrating factor 
R^(C i+) ,r) = 6Xp(r2) G C(r,exp(r 2 /2)), 

and the first integral 

r(+W+) r ^ = P;Pnr 2W {rK-B 2 P' n )h+A 
n ^ ' > (rP^- J B 2 P^)(exp(r 2 )+r^+ 2 (P n P^) 2 /i)' 



with 



J i= /^Q*. A = P' n ex P {r 2 )B 1 exp(J^\, 2* = r 2 - C (+) r - £ - 1, 
B 2 = r 2 + C«r - £ - 1, $ - eX P^)« 2r2 " 2 ^ " 2 ) P « " ^) 



P'n (exp (r 2 ) + r^P n PU^±dr 

Note that (C^ > r ) £ i* 1 , being .F an extension of C(r, adding the solutions of the 
Schrodinger equation = A\&. 

Example 3 (Coulomb). In this case we have V- = — 2 ^^ 1 ^ + 1 and we can take 

TjrH = r t+1 exp(-r), = - 1 G C(r). 

Note that due to Proposition 1 we have that DT is strong isogaloisian. Since 

DT(y.) = n=™ + 2 >-^ + i, 

we have that the potential V- is shape invariant. 
It is hold that 

'+ 1 \ 2 „ ] ( (1+ l x 2 



see [1, 3]. In particular, for A = 1 — ^ i+\+ n ) , A 7^ 0, we have 

\ t + 1 + n J ( ■ ' r P n t + l + n 



(+) _ e + i k e + i ex P yi+n+i 



4r— 2nr 



(n ' 2) " r P„ ^ + n + l o Wp2f ^(ii)' 

where P n are polynomials of degree n, which are related to the generalized Laguerre polyno- 
mials Ln \ for more details see [34]. 
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According to Lemma 1 for T(r) = £{£ + 1) — 2{£ + l)r+r 2 and N(r) = r 2 the vector field (15) 
becomes 

4- = (*+!)- Vf + 1> + (^) 2 ^ - -V*) ^ + r*£, 

and admits: 

The invariant curve 

/J-W-V) = -C W + — + P i - € C(cH >r ), 

r P n £+l+n ' 

with generalized cofactor 

^(C^r-) = -r 2 C« - (* + l)r - § r 2 + ^^r 2 e C(C<"\r). 
The exponential factor 

F (-) (C (-), r) = /+2Pnexp (z^:) eL2 , 

with generalized cofactor 

4 _) (C (_) , = r + {£ + l)r + ^r 2 - t+l r 2 G C(r). 
We note that the vector field X~ admits the generalized Darboux integrating factor 

/ 2(£+l)r \ 

4-»« W .0 = , eXP( ' Wn) ^ ^ C (cM, r ,exp (iitit 

7 P 2 r 2(f.+2)( (■(-) i g+1 i a t+l Y V U+l + n 

The first integral 



(« + n + l)r-^-'exp lg= 
^ ,(CH ' r) = I P„(( < + „ + l),.p. + (<+1 ) (< V + „ + i y )-^) + 1 I «P ( / * ) * * 



with 



A = £r(C H + 2) + nr(C (_) + 1) + r(C ( " } + 2)P n , 
T /exp(f^Ff) exp(^) 



2<+2p2 / 4r-2nr 



r 2l+2p2 



Note that 

(+) _ £ + 1 _ 31 + n + 3 (I + n + l)i* f 



r £ + n + l ((£ + n + l)i*-(^+l)P„) 

(^+l) 2 Pn 

{£ + n + !)((£ + n + - + l)P n ) ' 
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For suitability let assume fi = I + n + 1, thus we also have 

(+) exp (Jf) iU+*P nX AxK + B i p n + exp (4r) fi 2 (r - £ - 1) 



C(n,2) 



(exp (^) r2^+2p nX ( /U p^ _ (£ + l)P n ) + /Ue xp(4r)) 
with 

Ai = r^i*' + + - r(3^ + n + 3)), Bi = + l)(r(2£ + n + 2) - (£ + 



exp 



4r— 2nr 



After the Darboux transformation DT the vector field X n becomes 

w+) + + 2) 2(1 + 1) (+)2 \ a a 

The vector field admits the invariant curve 

/ ( + . (< ( + », r) ._ C ( + . + ^ + |__l±i_ £C(C ( + »,,. ) , 

with generalized cofactor 

ifW(C (+ V) = -C <+> - ^ - f + 6 C( C W.r). 

The generalized exponential factor 

i,< + )( C ( + ), r ) = / + ^ex P (=^) e ^, 
with generalized cofactor 

4+)(f(+ , r) ^ + |__l±l_ eC(r , 

Hence, the vector field admits the generalized Darboux integrating factor 

/ 2(l+l)r \ 



R (+) (£(+) r \ = 



2 : 



(£+l)r 
Tn 



and 

4 + )(C (+) ,r)Gc(c (+) ,r,exp(^ 
The first integral is given by relation 



Note that /^(C^j r ) £ ^\ being F an extension of C(r, C^) adding the solutions of the 
Schrodinger equation Bf 1 ^ = A\l/. 
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5 Final remarks 

This work is a first approach of applying Darboux transformation into polynomial vector fields 
that its associated foliation is of Riccati type. The strong isogaloisian property of the Darboux 
transformation guarantees that the transformed vector field is also of Riccati type. Additional- 
ly, we show that in the case of shape invariant potential is preserved the rational structure of 
the invariant objects such as the invariant curves, generalized cofactors, generalized Darboux 
integrating factors and Darboux first integrals. We remark that in this work we analyze the case 
of algebraically solvable potentials and in particular we consider potentials satisfying the shape 
invariance condition. A natural question arises: What happen in the case of quasi-algebraically 
solvable potentials? A more general study should be also done in the future for non rational 
potentials. 
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